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ABSTRACT 

'  This  report  investigates  the  relative  merits  of  a  delay-Doppler  imaging  radar  based  on  matched  filter  and 
Wigner-Ville  approaches.  Both  approaches  are  formally  equivalent:  the  relative  merits  of  each  method  are 
based  solely  on  implementation  issues.  Given  the  current  state  of  optical  delay-Doppler  »adar  and  signal 
processing  capabilities,  the  matched  filter  approach  provides  significant  advantages  over  a  Wigner-Ville- 
based  approach.  Additional  applications  of  the  Wigner-Ville  distribution  to  laser  radar  measurements  are 
discussed. 
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1.  INTRODUCTION 


Harper  Whitehouse  of  NOSC  proposed  an  alternative  delay-Doppler  imaging  algorithm  based  on  the 
Wigner-Ville  distribution  (WVDi.  Bill  Miceli  of  the  Office  of  Naval  Research  (ONR)  sponsored  an 
unconventional  imaging  program  at  Lincoln  Laboratory  to  evaluate  this  technique  relativ  e  to  more  traditional 
matched  filter-based  delay -Doppler  imaging  techniques.  This  report  presents  the  results  of  this  effort. 

The  primary  purpose  of  this  report  is  to  ev  aluate  the  relativ  e  merits  of  a  delay -Doppler  imaging  radai 
based  on  matched  filter  and  Wigner-Ville  approaches.  Both  approaches  are  formally  equivalent:  the  relativ  e 
merits  of  each  method  are  based  solely  on  implementation  issues.  Given  an  appropriate  transmit  waveform, 
the  output  of  a  practical  matched  filter  receiv  er  can  produce  a  delay  -Doppler  image  directly.  In  contrast,  a 
WVD  approach  requires  a  deconvolution  operation  to  obtain  the  desired  delay-Doppler  image.  For  nominal 
operating  conditions,  the  matched  filter  approach  can  process  data  over  a  relatively  small  delay-Doppler 
window  without  sacrificing  image  resolution.  This  can  significantly  reduce  the  computational  effort  required 
to  obtain  an  image.  To  maintain  image  resolution,  a  WVD  approach  must  process  data  over  a  time-frequency 
window  equal  to  the  time-frequency  extent  of  the  received  waveform.  This  requirement  becomes  prohibitive 
for  high  time-bandw  idth  product  waveforms.  Given  the  current  state  of  optical  range-Doppler  radar  and 
signal  processing  capabilities,  the  matched  filter  approach  is  preferable  to  a  Wigner-Ville-based  approach. 

While  the  Wigner-Ville  approach  does  not  appear  to  be  the  method  of  choice  for  delay-Doppler 
imaging,  it  may  be  useful  in  other  laser  radar  applications.  Specifically,  the  WVD  is  a  joint  time-frequency 
representation  of  a  signal  which  may  be  applied  to  the  analysis  of  time-varying  heterody  ne  or  autody  ne 
Doppler-resolved  measurements  of  vibrating  and/or  rotating  targets. 

This  report  is  organized  as  follows.  The  motivation,  definition,  basic  properties,  and  examples  of  the 
WVD  are  discussed  in  Section  2.  The  relation  between  the  Wigner-Ville  and  (symmetric)  ambiguity  function 
is  given  in  Section  3.  Section  4  is  a  discussion  of  the  application  of  both  the  WVD  and  ambiguity  function 
to  delay-Doppler  imaging.  Conclusions  and  alternative  applications  of  the  Wigner-Ville  distribution  to 
optical  delay-Doppler-resolved  laser  radar  measurements  are  presented  in  Sections  5  and  6.  respectively  . 
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2.  FORMULATION 


2.1  BACKGROUND 

Mixed  time-frequency  representations  attempt  to  simultaneously  analyze  the  energy  distribution  of 
nonstationary  signals  in  both  time  and  frequency.  These  techniques  may  therefore  be  interpreted  a'' 
generalizations  of  Fourier  analysis  techniques  which  determine  the  spectral  energy  density  of  stationary 
signals.  A  simple  mixed  time-frequency  representation  which  is  commonly  employed  is  the  short-term 
Fourier  transform  (STFT )  or  spectrogram.  In  this  technique  the  signai  is  assumed  stationary  over  some  short 
time  interval,  and  standard  Fourier  transform  techniques  (e.g..  FFT )  are  used  to  obtain  a  frequency  spectrum 
of  data  within  this  interval.  By  computing  spectra  over  successive  time  intervals,  a  time-frequency 
representation  of  the  signal  is  obtained.  In  the  context  of  laser  radar  measurements,  this  technique  has  been 
used  to  obtain  mixed  time-frequency  representations  of  heterodyne  and  autedyne  Doppler-resolved  mea¬ 
surements  of  rotating  targets.1'2  1  In  the  laser  radar  terminology,  this  mixed  time-frequency  representation 
has  been  referred  to  as  a  Doppler-time  intensity  (DTI)  plot.j  The  principal  shortcoming  of  STFT-based 
methods  is  the  trade-off  between  frequency  and  time  resolution.  In  addition,  the  STFT  does  not  reproduce 
the  energy  density  spectrum  and  the  instantaneous  power  of  a  signal.  A  class  of  time-frequency  represen¬ 
tations  w  hich  overcome  these  shortcomings  has  been  developed. ' °  One  such  representation,  the  Wigner- 
Ville  distribution  ( WVD).  was  first  introduced  in  the  area  of  quantum  mechanics  by  Wigner  in  1932  and  later 
into  the  field  of  communications  by  Ville  in  1948  (References  3  and  4). 

Joint  time-frequency  representations  of  a  waveform  have  properties  similar  to  those  of  joint  probability 
density  functions:  these  functions  must  be  positive  definite  and  satisfy  certain  zeroth-order  moment  or 
marginal  requirements.  Physically,  these  representations  estimate  the  time-frequency  energy  distribution  of 
nonstationary  signals  or  random  processes  and  must  therefore  satisfy  both  positivity  and  energy  conserv  ation 
constraints.6  For  a  signal  x(t)  with  Fourier  transform  X(w).  we  require  that  the  2-D  time-frequency 
representation  F(t.w)  have  the  following  properties: 


F(t,co)>0 

positivity 

(2.1) 

|  F(t,co)dco  =  |x(t)| 

instantaneous  power 

(2.2) 

J  F(t.to)dt  =|X((D)j” 

spectral  density 

(2.3) 

JJ  F(t.(i))do)dt  =  ||x|| 

total  signal  energy 

(2.4) 

Integrals  go  from  to  °°  unless  otherwise  noted.  While  the  WVD  satisfies  conditions  (2.2)  through  (2.4) 
and  provides  additional  useful  information  concerning  the  energy  distribution  of  the  function  in  time  and 


frequency,  it  does  not.  in  general,  satisfy  the  positiv  it y  condition  expressed  in  ( 2. 1 ).  Cohen  ct  ill.  has  e  pros  en 
the  existence  of  a  clas ,  of  positise  time-frequency  representations  which  satisfy  conditions  (2.1  i  through 
(2.4).  This  important  class  of  positise  time-frequency  distributions  is  not  discussed  here. 

2.2  DEFINITION  OF  THE  VVIONFR- VILFF  DISTRIBUTION 

The  cross  WVD  of  two  time  signals  x(t)  and  y(t)  has  the  following  bilinear  form: 

Wxv(t.(o)  =  Jx^t  +  ^)y^t-^)e"JtoTdT  .  (2.5) 

The  auto  WVD  of  the  time  function  x(t)  is 

Wx(*  m)  =  Wx x(t.co)  =  Jx^t  +  ^ jx*^t je_ja)Tdt  .  (2.6) 

There  is  a  functional  similarity  between  the  above  equations  and  the  (symmetric  lOrm  of  the)  cross  ambiguity 
and  auto  ambiguity  functions.  The  relation  between  the  ambiguity  function  and  WVD  w  ill  be  discussed  in 
Section  3. 

The  WVD  of  X(w)  and  Ytwi.  the  Fourier  transforms  of  x(t)  and  ytt).  respectively,  is 


wX.y(“-0  =  Jx^ai  +  |jY^a)-^  je  j;,d^ 


Cl)  —  — 


It  can  easilv  be  shown  that 


(2.7) 


Wx  Y(co.t)  =  Wx  y(t,w)  (2.8) 

In  words,  the  cross  WVD  of  the  spectra  of  two  signals  can  be  computed  direct) y  from  the  WVD  of  the 
respective  time  signals  simply  by  exchanging  the  frequency  and  time  variables. 

2.3  PROPERTIES  OF  THE  WKJNERMLLE  DISTRIBUTION 

A  partial  list  of  the  properties  of  the  WVD  is  provided  below.  Properties  and  proofs  may  be  found 
elsewhere  (see.  for  example.  References  7  and  8). 

2.3.1  Symmetry 


Wx  y(t.CD)  =  Wy  x(t.co)=>  Wx(t,w)  =  W*(t. ro) 
The  WVD  of  any  real  or  complex  function  is  real. 


(2.9) 
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2.3.2  Time  shifts 


A  time  shift  in  both  time  signals  results  in  an  identical  time  shift  in  the  cross  VVYD. 

Modulation 

Modulating  each  time  signal  by  exp  (  jilt  I  results  in  a  frequency  shift  of  the  resultant  cross  WYD  b\  11. 

2.3.4  WYD  for  Time-Limited  Signals 

If  the  time  signals  x(t)  and  ytti  are  restricted  to  a  finite  time  interval,  then  the  cross-WVD  distribution 
is  restricted  to  the  identical  time  interval  for  all  frequencies. 

2.3.5  WYD  for  Frequency-Limited  Signals 

If  the  spectra  of  signals  xttt  and  ytt  >  are  band-limited,  then  the  cross- VYYD  distribution  will  be  similarh 
band-limited. 

2.3.6  Nonlinearity 

The  WYD  is  a  bilinear  function  of  two  signals:  the  WYD  of  the  algebraic  sum  of  two  signals  is  not  the 
algebraic  sum  of  the  WYD  of  each  signal. 

Wx+).(t.O))  =  Wx(t.co)  + Wy(t.O))  +  2Re[wx  v(t,a')]  (2.10) 

The  third  term  in  Equation  (2.10)  represents  the  cross  or  interference  term  and  is  the  result  of  the  (second 
order)  nonlinear  form  of  the  WYD.  The  interference  term  may  pose  a  problem  in  the  analysis  ot 
multicomponent  signals.4  10  By  low -pass  filtering  the  VYYD.  one  may  reduce  the  interference  effects  h. 
sacrificing  resolution.  This  issue  is  similar  to  the  sidelobe-resolution  trade-off  in  standard  Fourier  analv  sis. 

2.3.7  Marginal  Relations  (Zeroth-Order  Moments) 

For  the  cross  WYD  we  have 

—  [  Wx  v(t,(0)dco  =  x(t)g*(t)  .  (2.11) 

2rt « 

Jwx  y(t.co)dt  =  X(w)Y*(co)  (2.12) 

which  reduce  to  the  following  expressions  for  the  auto  WYD 

■—  [wx(t.o))d(o  =  |x(t)|‘  ,  (2.13) 

2K  J 
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J  Wx(t.co)dt  =  (2.14) 

Finally  the  total  energy  in  Ml)  is  given  by  the  integral  of  the  auto  W'VD  over  the  entire  time-frequencv  plane 

^-ffwx(t,co)dtdco  =  J|x(t)fdt=||x|r  .  (2.15) 


2.3.8  First-Order  Moments 


The  marginal  relations  presented  in  Section  2.3.7  are  zeroth-order  moment  expressions.  Here,  we 
consider  first-order  moments  of  the  WVD  in  both  time  and  frequency.  The  W’VD  has  the  useful  proper)  \  that 
the  first-order  moments  in  time  and  frequency  are  the  signal  instantaneous  frequency  and  group  deiav. 
respectively . 

An  expression  for  U  (t).  the  average  frequency  of  the  W'VD  as  a  function  of  time,  is  given  by  the 
follow  ing  expression 


Q  (ts  JtoWx(t.(o)dto 
X  j  W;x(t.C0)d(0 
which  simplifies  tos 


n*(0  =  lm 


x'(t) 
x(t ) 


(2.16) 


(2.17) 


Here.  1m  |*|  is  the  imaginary -part  operation  and  the  prime  (')  denotes  differentiation.  This  expression  is 
identically  zero  for  real  valued  functions.  Let  xtt)  be  a  complex  analytic  signal  which  can  be  expressed  as 
follows 

x(t)  =  a(t)ei0(,)  .  (2.18) 


with  ad)  and  <f>tt  i  real  functions  of  time.  The  average  ‘'requeues  of  the  W'VD  for  such  a  function  is  given  by 

Ox(t)  =  0(t)  (2.19) 

The  availin'  frequency  of  the  It  l  I)  at  time  t  is  equal  u>  the  instantaneous  frequency  of  the  compiles  analvtie 
signal. 

An  expression  for  T  tu»).  the  average  time  of  the  W'VD  as  a  function  of  frequency,  is  given  by  the 
following  expression 


Tx(o>)  = 


}tW'x(t.(o)dt 

JWx(t.io)dt 


(2.20) 
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which  simplifies  as  follows* 


r 

Xx(co)  =  -Imj 


X'(to) 

X(co) 


(2.21) 


Here.  X(oj>  is  the  Fourier  transform  of  the  time  signal  \(t).  1m  |*J  takes  the  imaginary  part  of  the  operand, 
and  the  prime  denotes  differentiation.  This  expression  is  identically  zero  for  real  valued  functions.  Let  Xuo  i 
he  a  complex  analytic  signal  which  can  he  expressed  as  follows 

X(co)  =  a(o))e-'0(tl)l  .  (2.22) 

xxith  ado)  and  Bun)  .al  functions  of  frequency .  The  average  frequency  of  the  WVD  for  such  a  function  is 
given 

Tx(to)  =  -e(co)  (2.23) 

The  rime  average  of  the  U  l  l)  at  a  given  frequency  is  equal  to  the  Signal  group  dehr. . 


2.3.9  Second-Order  Moments 

Local  second-order  moments  in  time  and  frequency  (defined  in  the  usual  wax  i  are  a  measure  of  the 
spread  of  the  WVD  about  mean  time  or  frequency  values,  respectively .  As  tire  WVD  is  not  positive  definite, 
the  local  second-order  moments  of  the  W  VD  distribution  max  he  negative,  a  result  w  hich  is  inconsistent  x\  ith 
a  strict  probability  theory  interpretation.  Hoxvever.  these  moments  max  still  provide  useful  information:  a^ 
second-order  moments  tend  to  characterize  structural  features  in  the  WVD.  they  may  be  used  to  indicate  the 
structure  of  multicomponent  signals.1" 

2.4  ANALYTIC  SIGNAL  REQUIREMENT 

For  (complex)  analytic  signals,  the  first-order  moments  of  the  WVD  in  time  and  frequency  are  the 
instantaneous  frequency  and  group  delay  .  respectivjy.  Strictly  speaking,  these  properties  do  not  apply  to 
real  x  alueu  signals.  In  particular,  if  xve  let  the  signal  rtt )  represent  the  real  part  of  an  analy  tic  signal  xtt  i.  then 

Wr(t.(o)  =  -[Wx(t.(o)  +  W'x(t.-w)]  +  g(t.co)  .  (2.24) 

4 

where 

W  u.oj)  =  WVD  of  the  real  signal  m ) 

r 

W  (t.oj)  =  WVD  of  the  analvtic  signal  xtt; 

\ 

-dt.io)  =  oscillatory  fringe  pattern. 

Tne  WVD  of  a  real  signal,  as  depicted  in  Figure  1(b).  is  composed  of  three  terms,  the  desired  "image"  term 
W  (t. co).  a  conjugate  image  term  W  (t.-w).  and  a  residual  low-frequency  frmge  pattern  git.w.  Analysis  has 
shoxvn  that  the  artifacts  embodied  in  the  gu.w)  term  ha'  e  no  physical  meaning  (see.  for  example.  Refer¬ 
ence  9).  The  analytic  signal  requirement  must  be  take.i  into  consideration  when  implementing  the  W  VD. 
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WAVEFORM 


figure  /  Anulxtit  igmd  requirement .  The  time-frequency  <  huracter  of  a  linear  chirp  signal  is  given  in  (a)  Both  the  time  and 
frequent  \  axes  have  been  normalized  to  unit  tone  and  frequency.  respectively.  The  VV\  l)  of  a  real  signal  with  time-frequency 
i  harm  /eristics  given  in  *ai  is  shown  in  (hi.  The  U  \  I )  is  characterized  h\  the  image,  conjugate  image,  and  residual  fringe  pattern 
<  otnponents  The  U  \  I )  nf  the  amdvtu  linear  /  \f  chirp  is  shown  in  (e).  The  frequency  time  characteristics  of  the  signal  are  dearly 
reprodm  ctl  in  the  U  \  /) 
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2.5  MULTICOMPONENT  SIGNALS 


Laser  radar  measurements  of  complex  targets  are  generally  characterized  by  complicated  multicompo¬ 
nent  signals.  The  WVD  of  a  multicomponent  signal  is  complicated,  however,  by  two  factors:  ( 1 )  negative 
second-order  moments,  and  (2)  the  presence  of  interference  terms  (which  are  manifest  as  oscillatory  fringe 
patterns  in  the  WVD  —  see  Section  2.3.6).  Both  these  effects  may  he  reduced,  at  the  expense  of  resolution, 
by  low-pass  filtering  the  WVD  l>  Time-frequency  distributions  have  been  proposed  which  either  eliminate 
or  suppress  the  undesirable  properties  of  the  WVD.  >l0 

2.6  EXAMPLES 


In  this  section  we  demonstrate  the  basic  properties  of  the  WVD  by  example.  We  begin  by  considering 
the  WVD  of  a  simple  up-chirp  w  aveform,  as  shown  in  Figure  1(a).  The  WVD  of  the  up  chirp  is  show  n  in 
Figure  2(a).  The  time-frequency  characteristics  of  the  up  chirp  are  clearly  manifest  in  the  WVD.  The  first- 
order  moment  in  frequency,  which  is  the  instantaneous  frequency  estimate,  is  shown  in  Figure  2(b).  Note 
that  the  spread  or  variance  about  the  true  instantaneous  frequency  is  varying  but  hounded.  It  is  interesting 
to  note  that  the  local  variance  of  the  WVD  tends  to  increase  at  discontinuities  in  the  underlying  instantaneous 
frequency  of  the  signal. 


WIGNER-VILLE  DISTRIBUTION 
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FIRST-  AND  SECOND-ORDER 
MOMENTS  IN  FREQUENCY 
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Figure  J  t  he  H  )  l)  and  In,  a!  moments  <> /  a  single  up  i  lnt/>  I  lie  U  1 1)  of  il:e  tanalytu  i  up  <  Imp  is  sinus  n  in  iai  \arniahzed  time- 
frei/ueni  v  uses  nre  used  I  lie  In  si  and  set  on, I order  moments "/  the  U  \  It  nee  shim  n  in  /bias  the  solid  and  dashed  lines. respe,  lively, 
t  he  ineiin  represents  it  re,isomihle  estimate  < </  the  instantaneous  fiei/urney  ol  die  analun  up-,  Imp  waveform.  I  he  varianee  is 
reasonably  uniform  and  small,  in, In  tiling  that  the  "knife-edge'  ridge  defining  llie  ,  Imp  in  iai  is  well  defined. 
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The  WVDs  for  different  multicomponent  up-down  chirp  signals  are  shown  in  Figures  3  through  5.  The 
undesirable  oscillatory  fringe  pattern  in  these  figures  is  the  result  of  the  bilinear  form  of  the  WVP:  the  WVD 
of  the  composite  up-down  chirp  is  not  equal  to  the  sum  of  the  WVD  of  each  constituent  chirp.  Due  to  the 
oscillatory  nature  ol  the  cross-term  component,  this  term  may  be  reduced  by  low -pass  filtering  the  WVD. 

2.7  IMPLEMENTATION 

Although  computational  complexity  is  inherent  in  an  algorithm  which  generates  a  2-D  distribution  from 
1-Ddata.  implementation  issues  are  not  addressed  in  this  report.  Efficient  numerical  algorithms  may  be  found 
in  the  literature. ' !  !  Optical  implementation  of  the  W  VD  is  also  a  current  topic  of  investigation. ' '  1,1 
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3.  RELATION  BETWEEN  THE  WIGNER-VILLE  DISTRIBUTION 
AND  THE  AMBIGUITY  FUNCTION 


Both  the  WVD  and  the  ambiguity  function  tAF)  are  mixed  2-D  representations  of  a  I  -D  signal.  The 
utility  of  the  WVD  as  a  time -frequency  analysis  tool,  along  with  a  partial  list  of  applications  to  optical  radar, 
has  been  discussed.  In  the  context  of  optical  radar,  the  AF  is  a  useful  tool  to  evaluate  the  performance  of  a 
2-D  d-day-Doppler  matched  filter  receiver.  In  this  case,  the  magnitude  squared  of  the  ambiguity  function  may 
be  interpreted  as  the  point  spread  response  of  a  delay-Doppler  imaging  system. 

Fora  given  signal  x(t).  the  WVD  W  (t.wt  ann  me  ambiguity  lunction  A  (c.t>  are  defined  as  follows: 


Wx(t.{o)  =  Jx^t  +  |]x*(t-|je  JTUdT  .  (3.1) 

Ax(^,r)  =  Jx(^t  +  (3.2) 

These  equations  imply  that  the  AF  and  WVD  form  a  Fourier-transform  pair  1 

Ax(5.T)  =  FWx(t.co)  .  (3.3) 

Wx(t.(o)=F-'Ax(C.T)  .  (3.4) 

w  here  F  and  F  1  are  the  Fourier-transform  operations  defined  as  follows 

F  =  ~|J  e-j(^_UT)dtds  .  (3.5) 

F-1  =  JJej^I'£OT)dtd(o  (3.6) 


The  WignerA  ille  distribution  and  ambiguity  function  arc  two  different  representations  of  the  same  data.  A 
comparison  of  the  effect  of  certain  signal  operations  on  the  WVD  and  AF  are  listed  in  Table  1.  Note  that  time 
and  frequency  shifts  of  the  original  signal  are  preserved  in  the  WVD  representation,  but  are  mapped  into  a 
complex  oscillation  in  the  AF  representation.  Assuming  that  both  the  WVD  and  AF  have  similar  hardware 
(or  software)  implementations,  the  decision  to  use  one  method  over  another  depends  entirely  on  which 
method  provides  the  most  satisfactory  data  representation  fora  given  application. 

Figure  6  illustrates  the  ambiguity  diagram  and  WVD  for  an  up-down  chirp  w  aveform.  The  forw  ard 
Fourier  transform  of  the  product  x(t  +  t/2)  x  (t  -  t/2)  with  respect  to  delay  -  results  in  the  WVD.  w  hile  the 
forward  Fourier  tiansform  of  the  product  x(t  +  t/2)  x’(t  -  t/2)  w  ith  respect  to  time  produces  the  AF.  (The 
ambiguity  diagram  shown  is  the  magnitude  of  the  AF. )  The  2-D  Fourier  relation  is  clear:  the  inverse  Fourier 
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transform  of  the  WVD  with  respect  to  delay  followed  by  a  forward  Fourier  transform  with  respect  to  time 
produces  the  AF.  Note  that  time-frequency  characteristics  of  the  original  up-down  chirp  waveform  are 
clearly  evident  in  the  WVD.  The  fringe  patterns  evident  in  the  WVD  are  the  so-called  cross-terms  or 
interference  terms,  and  are  the  consequence  of  the  bilinear  form  of  the  WVD  kernel.  The  ambiguity  diagram 
of  the  up-down  chirp  is  characterized  by  a  central  peak  plus  relatively  high  sidelohes  w  hich  represent  the 
"ambiguity"  of  each  chirp.  Due  to  the  Fourier  relation  between  the  AF  and  WVD.  the  AF  becomes  more 
localized  in  the  delay  -Doppler  plane  and  the  corresponding  WVD  becomes  increasingly  distributed  through¬ 
out  the  time-frequency  plane.  This  is  an  important  result  in  the  context  of  delay -Doppler  imaging. 
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4.  DELAY-DOPPLER  IMAGING 


4.1  CONVENTIONAL  MATCHED  FILTER  APPROACH 

In  the  conventional  analysis  of'delay-Doppler  imaging,  the  output  of a  2-D  delay-Doppler  matched  filter 
represents  the  convolution  of  an  underlying  target  distribution  in  delay-Doppler  space,  denoted  by  itiCt) 
convolved  with  the  magnitude  squared  of  the  ambiguity  function  or  point  spread  response  (PSR )  of  the  filter. 
Unlike  conventional  spatial  imaging  systems  whose  PSR  depends  on  a  spatial  aperture,  the  PSR  of  the  delay  - 
Doppler  imaging  system  is  determined  by  the  time-frequency  properties  of  the  transmitted  waveform.  If.  for 
example,  one  could  generate  a  thumbtack-like  PSR.  the  output  of  the  2-D  matched  filter  would  represent  the 
2-D  delay-Doppler  image  of  the  target.  The  relation  between  the  output  of  a  2-D  matched  filter  and  the 
underlying  delay-Doppler  distribution  of  the  target  is  given  by  the  following  expression: Jl 

(|Ax.y(C- 4')  =  JJ°(C-T')|A*[C ~ (C-SoM - (T  - To) )J-dTdC  .  (4.1) 

where 

x(t)  =  transmitted  waveform 

y(t)  =  received  waveform 

a( £.7)  =  targe:  cross-section  distribution  in  delay  and  Doppler 

A  (£.t)  =  auto  ambiguity  function  of  transmitted  waveform 

A  (£.7)  =  cross  ambiguity  function  of  received  and  transmitted  w  aveforms 

t  =  filter  mismatch  in  delay  (tracking  error) 

£  =  filter  mismatch  in  Doppler  (tracking  error). 

Here,  the  operator  <>  indicates  an  ensemble  average  (e.g..  an  average  over  different  speckle  realizations). 
The  basic  assumptions  used  to  derive  Equation  (4.1 )  may  be  found  in  Rihac/.ek.20  The  output  of  the  2-D 
matched  filter  is  the  cross  ambiguity  function  between  the  transmitted  and  received  waveforms.  As  the  filter 
output  <IA  ^(£.7)l‘>  represents  the  convolution  between  the  desired  target  cross-section  distribution  (Kl.ti 
with  the  magnitude  squared  auto  ambiguity  function  IA(c.  7)1'  of  the  transmitted  waveform.  IA  (£.  7)1'  may 
be  interpreted  as  the  PSR  of  a  delay-Doppler  imaging  system. 

Given  flexibility  in  w  aveform  generation  and  reception,  one  could  select  a  waveform  w  ith  a  thumbtack¬ 
like  ambiguity  function  resulting  in  a  matched  filter  output  w  hich  approximates  the  underlying  target  cross- 
section  distribution;  there  is  no  need  to  invert  Equation  (4.1 )  for  rr(c.7). 

4.2  WIGNER-VTLLE  APPROACH 

Harper  Whitehouse  derived  a  delay-Doppler  imaging  technique  based  on  the  WVD."1  The  starting 
point  is  Equation  (4.1 ).  Substituting  the  following  two  relations14"  2' 

f|ax  v(^,t)!~  =  A*(i;.t)  Ay(^.t)  (4.2) 
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and 


F~*Ax(£.t)  =  Wx(t.a)) 
into  Equation  (4. 1 )  gives  the  following  result 


(4.3) 


(|wy(t.(o)^  =  Jja(t'.o)')Wx(t-t',(o-a)')dco’dt'  .  (4.4) 

This  expression  is  the  basis  for  delay  -Doppler  imaging  with  the  WVD.  It  states  that  the  expectation  of  the 
WVD  of  the  received  w  aveform  is  equal  to  the  convolution  of  the  underlying  delay-Doppler  target  cross 
section  with  the  WVD  of  the  transmitted  signal.  Note  that  ( 1 )  a  cross-WVD  computation  is  not  required.  (2 ) 
the  WVD  of  the  transmit  signal  may  be  interpreted  as  the  PSR  of  the  WVD-based  imaging  sy  stem,  and  (3) 
the  matched  filter  and  Wigner-Ville  imaging  systems  [as  defined  by  Equations  (4. 1 )  and  (4.4).  respectively] 
are  different  representations  of  the  same  underlying  process;  from  an  information  theoretic  point  of  view,  the 
methods  are  identical. 

4.3  COMPARISON  OF  THE  MATCHED  FILTER  AND  WIGNER-VILLE  APPROACHES 

As  previously  stated,  the  matched  filter  and  Wigner-Ville  approaches  are  equivalent  from  a  theoretical 
view  point.  Imaging  w ith  the  two  approaches  is.  however,  very  different.  The  comparison  of  the  matched 
filter  and  Wigner-Ville-based  delay-Doppler  imaging  systems  will  begin  w  ith  an  evaluation  of  the  point 
spread  response  (PSR).  The  PSR  is  a  fundamental  property  of  an  imaging  system.  The  computational 
requirements  for  each  system  will  then  be  discussed.  Bistatic  issues  will  be  addressed.  As  will  become 
evident.  S/'N  issues  are  secondary  and  w  ill  not  be  discussed  here.  A  preliminary  S/N  analy  sis  of  the  WVD 
may  be  found  in  References  23  and  24. 

4.3.1  Point  Spread  Response 

Thumbtack-like  ambiguity  responses  for  the  matched  filter  receiver  were  discussed  briefly  in  Sec¬ 
tion  4.1 .  In  the  case  of  waveforms  which  have  thumbtack-like  ambiguity  functions  (in  the  region  of  interest  i. 
delay  (range)  resolution  is  inversely  proportional  to  the  frequency  bandwidth,  and  Doppler  resolution  is 
inversely  proportional  to  the  temporal  duration  of  the  w  aveform.  Increasing  the  time-bandw  idth  product  of 
these  waveforms  allows  one  to  extract  an  increasing  amount  of  information  from  the  target.  The  time- 
bandwidth  product  of  a  signal  may  be  increased  by  phase  or  frequency-modulation  techniques.  These 
techniques  generally  lead  to  simple  implementations  relative  to  those  obtained  using  amplitude-modulation 
techniques.  The  main  point  is  that  it  is  possible  to  simultaneously  obtain  high  resolution  in  both  delay  and 
Doppler  by  choosing  a  waveform  with  a  thumbtack-like  ambiguity  function. 

The  WVD-based  receiver  has  significantly  different  behavior.  As  shown  in  Equation  (4.4).  the  PSR  of 
a  WVD-based  receiver  is  given  by  the  WVD  of  the  transmit  waveform.  A  thumbtack-like  WVD  requires  the 
transmit  signal  to  be  ( 1 )  nonzero  over  a  time  interval  small  compared  w  ith  the  target  extent  in  time,  and  (2) 
nonzero  over  a  frequency  interval  small  compared  with  the  frequency  extent  of  the  target.  These  two 
conditions  are  generally  mutually  exclusive;  a  thumbtack-like  WVD  may  be  unrealizable  for  many  targets 
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of  interest.  To  obtain  a  high-resolution  delay -Doppler  image  directly,  the  WVD  approach  requires 
narrowband  short-duration  waveforms,  while  the  matched  filter  approach  requires  wideband  long-duration 
waveforms.  As  wideband  long-duration  waveforms  are  physically  realizable,  the  matched  filter  approach 
has  a  significant  advantage  over  the  WVD  approach.  This  property  is  consistent  with  the  Fourier-transform 
relation  between  the  two  methods:  an  AF  which  is  highly  localized  transforms  into  a  WVD  which  is  highlv 
distributed.  This  is  demonstrated  for  an  up-down  chirp  in  Figure  b. 

4.3.2  Computational  Requirements 

We  begin  by  assuming  the  existence  of  waveforms  w  ith  thumbtack-like  ambiguity  functions.  While 
the  computations  involved  in  computing  the  full  matched  filter  output  and  WVD  output  are  identical,  the 
output  of  the  matched  filter  gives  the  desired  image  directly  whereas  the  WVD-based  method  requires  a 
deconvolution  to  obtain  the  image  o(c.t).  This  is  generally  an  ill-posed  and  computational!)  intensive 
operation. 

If  we  further  assume  that  both  the  delay-Doppler  tracking  errors  and  the  target  cross-section  delay  - 
Doppler  distribution  are  both  small  compared  with  the  delay-Doppler  extents  of  the  waveform,  then  the 
computational  requirements  of  the  matched  filter  approach  can  be  reduced  significantly  by  limiting  data 
processing  to  a  relatively  small  delay  -Doppler  w  indow  .  This  can  be  done  w  ithout  loss  of  image  resolution. 
From  a  practical  standpoint,  the  time-bandw  idth  capacity  of  the  signal  processor  may  be  reduced  w  ithout 
sacrificing  image  resolution.  In  order  to  maintain  image  resolution  using  the  W'VD  approach,  data  must  be 
processed  over  time-frequency  extents  defined  by  the  received  waveform.  As  a  waveform  with  a  time- 
bandwidth  product  equal  to  TB  has  a  2-D  W'VD  w  ith  (TBr  independent  resolution  cells,  large  TB  product 
w  aveforms  could  prohibit  real-time  computation  of  the  WVD. 

4.3.3  Bistatic  Considerations 

It  has  been  suggested  that  the  WVD  approach  may  have  some  advantages  in  a  bistatic  situation  as  the 
W'VD  receiver  produces  the  auto  WVD  of  the  received  waveform:  the  image  is  obtained  after  deconvolv  ing 
the  auto  W\  D  of  the  received  waveform  with  the  (stored)  auto  WVD  of  the  transmit  waveform.  In  this  case, 
a  reference  waveform  generator  is  not  required.  This  reasoning  is  somewhat  misleading:  the  transmit 
waveform  could  be  "stored"  in  a  programmable  matched  filter  or.  if  the  receive  and  transmit  waveforms  could 
be  digitized,  a  completely  flexible  matched  filter  could  be  implemented  digitally  .  While  these  issues  are 
technology  driven,  it  does  not  appear  that  the  WVD  approach  would  provide  significant  advantages  over  a 
matched  filter  approach  in  the  case  of  bistatic  measurements. 
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5.  CONCLUSIONS 


The  primary  purpose  of  this  report  is  to  evaluate  the  relative  merits  of  a  delay -Doppler  imagine  radar 
based  on  matched  filter  and  Wigner-Ville  approaches.  Both  approaches  are  formally  equivalent:  the  relath  e 
merits  of  each  method  are  based  solely  on  implementation  issues.  Given  an  appropriate  transmit  waveform, 
the  output  of  a  practical  matched  filter  receiver  can  produce  a  delay-Doppler  image  directly .  In  contrast,  a 
WVD  approach  requires  a  deconvolution  operation  to  obtain  the  desired  delay-Doppler  image.  For  nominal 
operating  conditions,  the  matched  filter  approach  can  process  data  over  a  relatively  small  delay-Doppler 
window  w  ithout  sacrificing  image  resolution.  This  can  significantly  reduce  the  computational  effort  required 
to  obtain  an  image.  To  maintain  image  resolution,  a  WVD  approach  must  process  data  over  a  time-frequency 
window'  equal  to  the  time- frequency  extent  of  the  received  waveform.  This  requirement  becomes  prohibitiv  e 
for  high  time-bandw  idth  product  waveforms.  Given  the  current  state  of  optical  range-Doppler  radar  anti 
signal  processing  capabilities,  the  matched  filter  approach  is  preferable  to  a  Wigner-Ville-based  approach. 
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6.  FUTURE  WORK 


In  optical  radar  applications.  time-frequency  representations  have  been  used  to  process  both  heterody  ne 
and  autodyne  Doppler-resolved  cross-section  measurements  of  rotating  and/or  vibrating  targets.  In 
particular,  the  short-term  Fourier  transform  (STFT i  has  been  used  to  analy  /e  the  time-frequency  distribution 
of  both  autodyne  and  heterodyne  laser  radar  measurements  of  rotating  targets.  Standard  ('requeues 
discrimination  techniques  hav  e  been  applied  to  the  analy  sis  of  laser  radar  Doppler-resolved  data  of  vibrating 
targets.  An  investigation  into  the  relative  merits  of  standard  time-frequency  analysis  techniques  (e.g..  STFT. 
frequency  discriminator,  etc.)  vs  signal  analysis  using  time-frequency  distribution  techniques  (e.g..  the 
WVD)  may  prove  useful  in  these  applications. 
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